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1. INTRODUCTION 

The main issues in robot control problems include accurate reference tracking, load disturbance re- 
jection and robustness to parameter variations. When all these relatively important requirements are taken 
into account, the high-precision robot control problem presents some challenges due to the coupled nonlinear 
dynamics of robotic systems. Various control schemes that may be divided into model-based and model-free 
control categories have been designed and implemented for robot manipulators to tackle with these challenges. 
The model-based controller designs include sliding mode control (SMC) [I], [2], computed torque control [3], 
robust control [4], adaptive control [5], model predictive control [6] and backstepping control [7]. The model 
uncertainties stemming from frictions, parametric uncertainties and load disturbances are the inherent issues 
of the model-based controllers. To eliminate the effects of these uncertainties, robust, adaptive and machine 
learning techniques have been used in the model-based control designs [8]-[10]. Specifically, many forms of 
SMC approaches have been designed for robot manipulators including traditional SMC [LI], integral SMC 
[12], fuzzy SMC [13], backstepping SMC [7], [14], terminal SMC [15], [16], second-order SMCs [17], non- 
singular terminal SMC [18], fast terminal SMC [19], time-delay estimation based SMC and adaptive SMC 
[21], [22], but these approaches require an accurate robot model to assure a good control response. 

In recent years, some model-free control approaches have also been designed for robot manipula- 
tors. Since, in general, a model (or model accuracy) is not considered in the control designs, the model-free 
control methods offer some advantageous alternative approaches for robot control problems. The PD and PD 
plus gravity controllers are the most common and highly effective model-free or minimum model requiring 
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controllers for the industrial robotic manipulators [23], [24]. Machine learning-based controllers including 
neural networks and fuzzy logic approaches are also designed in this context [25]. A data-driven observer- 
based model-free adaptive discrete-time terminal SMC [26], an extended state observer-based controller [27], 
the PID control 28], and an intelligent proportional-derivative SMC have also been designed as 
model-free controllers. The designed model-free controllers mostly consist of a controller plus estimated robot 
dynamics, i.e., the alternative model-free controllers have virtually the same structure as the PD plus gravity 
controller, but they have more complicated control design procedures and their control performances are not 
significantly better than the PD plus gravity controller. 

In this study, a model-free integral SMC is proposed for robust control of robot manipulators without 
any need for the system model. The inherent robustness of SMC can properly cope with the highly coupled non- 
linear dynamics of robot manipulators. Many SMC forms have been designed for robot manipulators, but these 
approaches require a robot model to assure a good response or otherwise can exhibit poor tracking performance. 
In this work, a three-term model-free continuous integral SMC is designed in order to assure high tracking 
accuracy, fast convergence, and strong robustness against disturbances. The input—output information of the 
robot without using any model parameters is utilized in the proposed controller design. The numerical results 
in both joint-space and task-space implementations are provided for the PUMA-560 robotic arm to illustrate 
the performance of the proposed controller. The results of the method are compared with the industrial PD plus 
gravity controller. The proposed control method is able to provide a short settling time, non-overshoot response, 
zero steady-state error, and robustness under the existence of load disturbances and parameter variations. 

The study is organized as follows. A model-free integral SMC design for both joint-space and task- 
space controls is presented in sections 2 and 3, respectively. Section 4 presents an application of the control 
method to the PUMA 560 robot, and section 5 concludes the paper. 


2. JOINT-SPACE ROBOT CONTROL WITH A MODEL-FREE INTEGRAL SMC 
The closed-form dynamic equations of an n-DOF robotic manipulator can be given by (1): 


M(qg)é+C(a,aa+G(q) =7T (1) 


where the joint parameters are defined with n-dimensional joint position vector q(t), joint velocity vector q(t) 
and joint acceleration vector g(t). The actuation torques for each joint are given by 7 € R”. The positive 
definite and symmetric matrix M(q) € R”*” is the inertia matrix, C(q, q) € R”*” is the Coriolis/centripetal 
term, and G(q) € R” is the vector of gravitational forces. The inertia matrix is bounded with mI < M(q) < 
My where m; and m,, denote positive constants and J is the identity matrix. The robotic manipulators satisfy 
the skew-symmetric relationship, i.e., 


q' (M(q) — 2C(q,4))q = 0 (2) 


For all g € R”, the Coriolis/centripetal term is bounded by a quadratic function as (3): 


IIC(4, aall < collall” (3) 


where co is a known scalar function (or a positive constant). The gravity vector is also bounded with ||G(q)|| < 
go for a finite constant go > 0. 

The control purpose of the robotic manipulator is trajectory tracking control (or reference regulation 
as its special case). Designing joint-space controllers with smooth joint-space trajectories is the most common 
approach. Consider a tracking error vector with a position control objective as (4): 


where qq is a reference signal vector. The control problem is to find a suitable actuation torque 7 such that the 
position error goes to zero, e > Nas t — oo. The control system must achieve this control objective with a good 
transient response to reference changes, and it should reject load disturbances and tackle measurement noise. 
The common industrial controllers include the PD control or PD plus gravity control methods for providing 
desired control goals for robot manipulators. Hence, the PD plus gravity controller will first be revisited, then 
a model-free controller will be introduced, and finally, the results will be provided comparatively. 
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2.1. Revisiting the industrial PD plus gravity controller 

Reference regulation of robotic manipulators is commonly accomplished by independent joint PD 
control law with a gravity compensation term [23]. It is well-known that a full compensation of gravity vector 
is not fully possible for degree-of-freedom (n-DoF) rigid manipulators, while the control performance can 
significantly be improved with the gravity term. The PD plus gravity controller has the form of (5). 


7 =—Kye— Kae + G(q) (5) 


where the control matrices K, € R"*” and Kg € R”*” are positive-definite and diagonal matrices and 


G € R” is an approximate value of the gravity vector G (q). Substituting (5) into {i}, the closed-loop system 
can be written as (6): 


M(q)G+ C(4,4)4 + Gq) + Kae + Kye = 0 (6) 


where G = G — G. If the reference is a constant, then the PD plus gravity controller is quite effective. For 
stability analysis, consider a positive-definite Lyapunov function as (7): 


1 
2 


The time-derivative of (7) is negative definite for a constant reference signal gq = 0 and the skew-symmetric 
robot property g7 (IM — 2C)q = 0, namely, 


1, . 
V = 54° M(a)at 5e" Kye %) 


V < -Amin(Ka)llall? + Glial (8) 


where \;nin is the smallest diagonal element of Ky. Thus, the stability is assured with V <0 i. (Ka)|lq|| > 
||G/|. The steady-state error approaches to zero, but not becoming zero because the equilibrium point of the 
closed loop system re f6 is found as (9): 


G(de) + Kpee = 0 > €e = G(ae)/Kp (9) 


where qe and €, denote equilibrium points. Since ||G|| < go for a constant go > 0, the error is bounded. 
In addition, in practical implementations, the term G(q) is replaced with G(qa) (known as PD plus desired 
gravity compensation), where qq is the reference signal. Therefore, the PD plus gravity controller requires 
large proportional control gains to drive the error to the close vicinity of zero, which might excite measurement 
noise and high-frequency system dynamics [23]. 


2.2. A model-free integral SMC design 

The SMC provides robust control solutions, but the reaching and sliding phases may require a long 
time and high control effort due to the conservative bounds of the model uncertainty. To ensure high control 
performance, fast convergence, and robustness against disturbances and robot manipulator variations, a three- 
term model-free integral SMC can be defined by (10): 


T = —Kos — Kj|s|“sgn(s) — 70 
79 = Ko|s|%sgn/(s) (10) 
s=e+ Ke 
where tracking error e is defined as e = q — qq for a reference signal gg, control gain matrices are de- 


fined with positive diagonal components as Ky = diag{ko1,...,kon}, Ki = diag{ki,..., kin}, Ko = 
diag{k21,..., kan} > 0, K = diag{ki,...,kn} > 0, and 0 < a < 1. The function, |s|%sgn/(s), is defined in 
the component-wise form as (11): 


|s|“sgn(s) = [|s1|*sgn(s1),---,|Snl*sgn(sn)]” (11) 


where sgn(s;) = s;/|s;|. To facilitate the operations, a virtual reference vector q, is defined as (12): 


dr =Q-8 (12) 
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Substituting the control law into the robot equation (i), the closed-loop system dynamics are obtained as 
(13): 


Mé =r — MG, — CG, — G(q) — Cs (13) 


Consider a positive definite Lyapunov function for stability analysis of the closed-loop system as (14): 


ve so" M(q)s (14) 


The time-derivative of can be written as (15): 


: 1 : 
V=s'Mét+ 5° Ms 


1 ; 
= s"lr— Md, — Cd, — G] s'Cs+58'Ms 


n 1 

=—s'Kos—s™A So kulsilet* — 
i=l 

< —S fkoslsi| + kys|si|* — |Aal]| sal 


i=l 


where s™(M — 2C)s = 0 due to the skew-symmetric property, and A = Mq, + Cq, + G + 79. The stability 
is achieved with V < 0 if (ko;|5;| + k1|5;|%) > |A,|. For a = 0.5, a solution of this condition yields (16): 

1 
ORR, 
where |A;| < do. Thus, the control signal drives the robot states into the bounded set in finite time. Since 
V <0, the system is stable, and the trajectory reaches the sliding surface in finite time and stays in a very close 
neighborhood of it. The bounded set can be designed to be very small with large control gain choices, i.e., the 
robustness and practical stability of the controller are assured when k1; > do and ko; > 0. Once the robot 
trajectory reaches this bounded set, then all the state variables will virtually have their steady-state values. For 


a constant reference signal qa, a unique steady-state solution of the closed-loop system (13) is found by setting 
all derivatives to zero as (17): 


|s;| > [k2, + Qkoido — k1i(k?, + 4kosdo)!/?] (16) 


G(qde) = -K, | |s-|°sgn(s_)dt (17) 
0 


where qe and s- denote equilibrium points. This implies that the term A in will approach zero, dg — 0 
with ¢ — oo, and thus, the bounded set will go to zero |s;| — 0. Due to the fact that the differentiable 
function V(t) has a finite limit as t + oo, and V(t) is uniformly continuous, V(t) + 0 as t — oo according 
to Barbalat’s lemma [31]. This means that the error functions e and é tend to zero as time goes to infinity. 

There are five control parameters to be determined. The features of these control parameters are as 
follows. The positive diagonal control gain matrices Ko and KK; must be selected large enough to guarantee 
system stability. The integral term of the controller acts as an adaptive controller as seen in equation (17). and 
is used to compensate for the static system parameters and disturbances. Hence, the disturbance compensation 
speed of the controller is determined from the value of the gain matrix K2, which can be selected as ko; > go 
(upper bound of gravity vector). The control parameters K determine the speed of the transient response of 
the control system. The parameter a is selected as a = 0.5 for good robustness and noise rejection features. 
Finally, the controller may encounter the chattering phenomenon, so the discontinuous switching function, 
sgn(-), may be replaced with its smooth approximates (e.g., sigmoidal or saturation functions) in practical 
applications with guaranteed uniform ultimate boundedness of errors. 


2.3. Task-space robot control with model-free integral SMC 
The joint-space control method describes robot movements using the torque and angular positions and 
velocities that are necessary to complete the task. While the joint-space control is designed by default, it needs 
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inverse kinematics to calculate the desired joint-space reference trajectories for completing a task, which pos- 
sesses some difficulties. An alternative approach is the task-space (or Cartesian-space) control, in which 
the forward kinematics and Jacobian matrix calculations are used to transform between joint-space variables 
(q1,--+-;@n) to task-space variables (21,..., 2m). The task-space control is commonly designed with the in- 
verse Jacobian or transpose Jacobian regulators to deal with a mechanical impedance of the robot end effector 
(32), (33). The advantage of the task-space control approach is to operate directly on the task-space variables, 
but the effects of singularities and redundancy cannot easily be managed, and it may become computationally 
demanding if velocities and accelerations are of concern [34]. 

The forward kinematics is used to find the pose of the end-effector. Given the joint angles gq € R”, 
the end-effector pose x € R™ is defined by (18). 


x= f(q), f: RB" R™ (18) 


The differential relationship between the joint displacements and the resulting end-effector motion is repre- 
sented by the Jacobian matrix as (19): 


z= J(q)d, J(a) = Of (q)/0¢ (19) 


where J(q) € R™*” is the analytical Jacobian matrix standing for the mapping from joint space to task space. 
The transpose of the Jacobian is also used to describe the torque and force relations by (20), 


7T=J(q)'F (20) 


where F' € R’ is the applied force at the end-effector in task space. The Jacobian transpose method eliminates 
the issues related to the Jacobian inversion and singularity problems, and thus it is used to obtain the robot 
model in the task space as (21), 


M(x) + C(a,4)¢ + G(x) =F (21) 


where M(x) = J-T M(q)J—, C(z,#) = J-TC(q,¢)J—! — J-TM(q)J-1 J J—! and G(x) = J-TG(q). 

It is assumed that the Jacobian matrix is nonsingular, but if it is singular then the system has re- 
dundancy for the given task in the task space. The task-space robot equation has the similar properties to 
joint-space robot equation, including the skew-symmetric property aT (M —2C)« = 0 and boundedness of the 
other robot equation terms [35]. Therefore, the joint-space controllers can also be designed in a similar form 
for the task-space representations. 

In task space a tracking error vector with a position control objective is defined as (22): 


{e=x-2a: jim e=o} (22) 


where xq is a reference signal vector. To achieve this tracking control goal for (2p, similar to a model-free 
integral SMC can be defined as (23): 


F = —Kos — Ki |s|“sgn(s) — fo 
fo = Ko|s|“sgn/(s) (23) 
s=e+Ke 


where e = x — &gq for a reference signal xg. Similar to the joint-space controller analysis, for a Lyapunov 
function, V = 0.58? M(zx)s, we can show that its time-derivative is negative definite as V < 0 under the 
controller (23). Namely, the error functions e and € tend to zero as time goes to infinity. 

The general task-space control and joint-space control schemes are illustrated in Figure 1. There are 
different pros and cons to trajectory planning and control in joint space or task space. For example, joint-space 
trajectories have smooth actuator motion and have no problem with singularities, while inverse kinematics are 
needed to generate joint-space references. On the other hand, task-space trajectories are easier to visualize and 
the motion in task-space is more predictable, but it is prone to singularities and the actuator motion is more 
difficult to verify. 
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Figure 1. Model-free control implementations in joint-space and task-space of robotic manipulators 


3. CONTROL IMPLEMENTATIONS TO PUMA 560 ROBOT 
The explicit dynamic model of the PUMA 560 robot is given by [36]: 


M(g)Gé+ B()adG +CladG+Gq=T-71 (24) 


where 7; is the load torque, B(q) is the matrix of Coriolis forces, C'(q) is the matrix of Centrifugal forces, and 
i # j fori,j = 1,...,6. The PUMA 560 robot matrix/vectors can be defined by M = [m,,], B = [bi;x], 
C = [ci], G = [Gi], whose components are given by [36], 

miu = Iml+h+13CCo+ I7+lo1+2h5)S'$23+ (ho +21e2)S'C23 +111 SC2+2C2[I5 $23 + 11223 + 116 $23] 
myz = IgS2 + IgCo3 + Ig9C2 + 113523 — ligC23 

my13 = IgCo3 + 11353 — I1gCo3 

m2 = Im2+ Ip + Ig + 2[15S3 + IeC2 + lis + Le $3] 

m3 = I5 $3 + Ig + Ty2C3 + T1683 + 2hhs 

m33 = Im3 + Ie + 2N5 

m35 = Ts + Liz 

maa = Im4 + Tha 

M55 = Im5 + Tz 

Me6 = Im6 + 123 

m1 = M12, 731 = 773, M32 = M923 

big = 2[—13S C24 U5 +lig) C2234 (I7+2l15 +121) SC23— 112 $293]+ (110 +2122) (1-253) +411 (1-282) 
bi13 = 2[(U5 + Li6)C2C23 + U7 + 2Ih5 + Io1)SC23 — Iy2C2S23 + I22(1 — 25$23)] + To(1 — 25993) 

bits = 2[-—SCo3 + Lys SC23 + I16C2C23 + I22CC 23] 

bi23 = 2[—IgS23 + I13C23 + tigS23] 

be1a = 114523 + [19523 + [20523 

bo23 = 2[—N2S3 + I5sC3 + IeC3] 

boas = 2[16C3 + 122] 

bo35 = 2[fi6C'3 + Ira] 

b314 = [2023 + traSo3 + [i9$23 


ba12 = —be1a4 
ba13 = —b314 
bais = —L20S23 — [17 S23 
bs14 = —b415 


C12 = I4Cz — IgSo3 — I9S2 + 113C23 + igS23 
C13 = 0.5b193 


C21 = —0.5b112 
€23 = 0.5b223 
c31 = —0.5b113 
C32 = —C23 

€51 = —0.5b115 
€52 = —0.5b295 


Gg = 91 Cz + g2S23 + 93S2 + gaC23 + 95523 
G3 = 92523 + gaC23 + 95523 
G5 = 95523 
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where S$; = sin(qi), C; = cos(q:), Ci; = cos(qi + qj), Sijk = sin(qi + g; + a), CC; = cos(q;)? and 
C'S; = cos(q;) sin(q;), and all other elements of the matrices are zero. The PUMA 560 robot arm constants 
and DH parameter values are provided in Tables 1 and 2, respectively. 


Table 1. Inertial (kg/ m?) and gravitational (Nm) constants 
No Inertial (kg/ m2) — Gravitational (Nm) constants No _ Inertial (kg/ m?) — Gravitational (Nm) constants 


1 ia 1.43 18 Tis 0.000431 
2 Ib 1.75 19 Tig 0.0003 

3 Is 1.38 20 Too —0.000202 
4 I 0.69 21 Ini —0.0001 
5 Ts 0.372 22 Too —0.000058 
6 Ig 333 23 Tog 0.00004 
7 Ir 298 24 Imi 1.14 

8 Is —0.134 25 Tes 4.71 

9 Io 0.0238 26 Im3 0.827 
10 Tio —0.0213 27 Ima 0.2 

ll Ly —0.0142 28 Ins 0.179 
12 Tho —0.011 29 Imé 0.193 

13 ls —0.00379 30 n —37.2 
14 la 0.00164 31 go 8.44 
15 Tis 0.00125 32 93 1.02 

16 lig 0.00124 33 4 0.249 

7 Tye 0.000642 34 95 —0.0282 


Table 2. PUMA 560 arm DH parameters 


Link 7 "di ay ay d; 
1 nm -90 0 0 
2 qa 0 ~=—0.4318 ~—-0.14909 
3 qa -—-«90-——*0.0203 0 
4 qa -90 0 0.43307 


The PUMA 560 robotic arm has six axes, but the last three axes form only a spherical wrist as seen in 
Figure 2. The control results are provided when the last three axes (q4, qs, q6) of the robotic arm are locked, 
and thus results corresponding to the first three axes (q1,q2,q3) are provided. The torque constraints of the 
PUMA 560 robot are taken into consideration in joint-space and task-space control designs, which are given 
by |7| < 97.6, |72| < 186.4, |73| < 89.4, |74] < 24.2, |75| < 20.1 and |7%| < 21.3 Nm. In numerical studies, 
the contributions from measurement noise and load disturbances are considered. It is assumed that the angular 
position and angular velocity measurements are exposed to normally distributed measurement noise with a 
mean of 4 = 0 and a standard deviation of o(q) = 0.172 deg (or +0.5 deg in magnitude) and o(q) = 0.5 deg/s 
(or +1.43 deg/s in magnitude). Step load disturbances with 7; = 10 Nm magnitude (around 10% of the first 
three joint torque limits) are added to system torque inputs at ¢ = 5 s for evaluating load disturbance rejection 
performances of the controllers. 


Figure 2. Schematic diagram of the PUMA 560 robotic manipulator 
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3.1. Joint-space control results 

Joint-space control results are provided in Figure 3 to Figure 8 for both PD plus gravity and the 
proposed model-free integral SMC controllers. The control gains are selected such that controllers provide 
non-overshoot responses and take into account the torque constraints of the PUMA 560 robotic manipulator. 
The PD plus gravity controller gain matrices are taken as K,, = 500 I3 and Kg = 100 J3 where [3 € Hees 78 
the identity matrix. The model-free integral SMC controller gain matrices are chosen as Ky = 60 I3, Ky = 50 
Is, Ko = 200 Is, K = 5 I3 and a = 0.5. It should be noted that the discontinuous function of the controller is 
approximated as sgn(s;) ~ sat(s;/5) in the implementations to avoid possible chattering in the torque control 
inputs. The control parameters are selected such that both controllers can have very close performances for 
comparison targets. 

Performances of the controllers for constant references, gq = (7/4,7/2,—7/4), are illustrated in 
Figures 3 to 5. Figure 3 shows that both controllers have non-overshoot transient responses and reach the 
steady-state conditions in less than 0.5 s with zero steady-state error under no-load conditions. Figure 4 shows 
the torque control signals of the controllers, in which both controllers produce almost the same torque inputs. 
Load disturbance rejection performances of the controllers are assessed for a step load 7, = 10 Nm (less 
than 10% of the joint torque limits) added to system torque inputs at t = 5 s as seen in Figure 4. Under the 
load disturbance, Figure 5 shows that the model-free integral SMC controller quickly eliminates the effects of 
the load disturbance and provides zero steady-state errors. However, the PD plus gravity controller provides 
an erroneous performance under load disturbance. Hence, the proposed model-free integral SMC controller 
provides robustness to the robotic manipulator control against load disturbances. 

Figure 6 shows the disturbance estimation performance of the controller. The integral term 7) of 
the controller acts as an adaptive controller and compensates for the static system terms and disturbances. 
Under no-load conditions, it is seen in the figure that the integral control components approach to the gravity 
components, 791 + G1, To2 + G2, and 73 — G3. When a step load 7; = 10 Nm is added to system torque 
inputs at f = 5 s, the integral control components approach to the total values of static terms, i.e., 71 — 
Gi +71, T2 4 Go +Ti2, and 73 + G3 +73. The gain matrix Ky» determines the disturbance compensation 
speed of the controller, i.e., the larger the gain, the faster convergence is observed, but an oscillation in the robot 
trajectory can be seen with large integral control gain values. Hence, there is a trade-off between controller 
performance and disturbance rejection response speed. 


PD plus gravity controller 
2 ie - 
Q 
q 8 0.5 
= = 0 YN 
= > 
S £05 . 
7 8 | 43 
1 . n L 1 1 1 1 f 
0 2 4 6 8 10 0 2 4 6 8 10 


tis) 


Model-free integral SMC controller 
T T 4 | 7 a r 


0.5 T)=10 Nm 


Joint Angles (rad) 


-0.5 


Joint Velocities (rad/s) 
oO 


“0 2 4 6 8 10 “0 2 4 6 8 10 
t(s) t(s) 
Figure 3. Constant reference tracking control performances for joint positions and velocities under 7; = 10 
Nm load torque 
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PD plus gravity controller Model-free integral SMC controller 
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ae T,=10 Nm 
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Figure 4. Control torque signals for each joint for constant references 
PD plus gravity controller Model-free integral SMC controller 
0.2F 0.2F 
e 
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Figure 5. Constant reference tracking errors under 7; = 10 Nm load torques 
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Figure 6. Disturbance estimation with an integral term of the model-free integral SMC controller 


Figures 7 to 9 show the performances of the controllers for sinusoidal reference signals. Figure 7 
displays non-overshoot transient responses and less than 0.5 s settling time with zero steady-state error in joint 
positions and velocities for both controllers under no-load disturbances. When a step load 7; = 10 Nm is added 
to system torque inputs at t = 5 s as seen in Figure 8, the PD plus gravity controller provides a constant steady- 
state error. On the other hand, the proposed model-free integral SMC controller quickly recovers the effects of 
load disturbances. The load disturbance rejection performances of the controllers can easily be seen in Figure 
9. The phase-plane portraits (x vs y) are plotted when a step load 7; = 10 Nm is added to system torque inputs 
at t = 0s as seen in Figure 9. It is clear that there is an offset error for the PD plus gravity controller. Similar 
to the constant reference case, the model-free integral SMC controller rejects the load disturbance and provides 
zero steady-state error and robustness to the robotic manipulator. 
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Figure 7. Sinusoidal reference tracking control performances for joint positions and velocities under 7; = 10 
Nm load torque 
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Figure 8. Control torque signals for each joint for sinusoidal references 
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Figure 9. Sinusoidal reference tracking errors under 7; = 10 Nm load torques 
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3.2. Task-space control results 

For task-space control design, the forward kinematics and its Jacobian matrix are needed as seen in 
Figure 1. Forward kinematics of the PUMA 560 robot obtained from the DH parameters given in Table 2 when 
the last three axes (qu, qs, qo) are locked is given by (25): 


CGO ein HS) 08, 
y = Sy(agCz + a3C23 + d4$23) + d2C1 (25) 
Z = A252 + a3S93 — d4C23 
where S; is used for sin(q,), C1 is used for cos(qi), S23 is used for sin(q2+q3), and C23 is used for cos(q2+4q3). 
The robot parameters are dz = 0.1491, d4 = 0.4331, ag = 0.4318 and a3 = 0.0203. The Jacobian matrix of 
the robot is obtained as J = [j;;] € R°*? with (26). 
jus = —S1(a2C2z + a3C23 + d4S23) — d2C1 
jiz = C1(—a2$2 — a3S23 + d4Co3) 
jiz = C1(—a3$23 + daCo3) 
jar = C1(a2C2 + a3C23 + d4S23) — d2Sy (26) 
j2z = S1(—a2S2 — a3$23 + d4C23) 
j23 = S1(—a35$23 + daC23) 
js1 =0,  js2 = d4S23 + a3C23 + a2C2 
Jj33 = d4S23 + a3C23 


Task-space control results for the PUMA 560 robotic manipulator are shown in Figures 10 to 12 for both PD 
plus gravity and the proposed model-free ISMC controllers. The PD plus gravity controller gain matrices 
are taken as K, = 5000/3 and Ka = 1000/3 where Jz € R3** is the identity matrix. The model-free 
ISMC controller gain matrices are chosen as Ky = 6001/3, K, = 500/3, Kg = 2000/3, K = 5/3 and 
a = 0.5. To eliminate possible chattering in the torque control inputs, the signum function is approximated 
as sgn(s;) * sat(s;/5) in the control implementations. Again, as in joint-space control designs, the control 
parameters are selected such that both controllers can have very close performances for comparison goals. In 
addition, non-overshoot responses for both controllers are targeted. 
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Figure 10. Task-space tracking control performances for joint positions and velocities under 7; = 10 Nm load 
torque 
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Task-space references are taken as xq(t) = 0.25+0.1 sin(¢), ya(t) = 0.25+0.1 cos(t) and za(t) = 0. 
Figure 10 displays task-space tracking control performances for joint positions and velocities under 7; = 10 
Nm load torque applied at t = 5s. It is obvious that both controllers have non-overshoot transient responses and 
reach the steady-state conditions in less than 0.5 s for both position and velocity reference signals. Figure 11 
shows that both controllers produce almost the same torque control signals. Under the load torque disturbances, 
Figure 12 clearly shows that the model-free integral SMC controller quickly rejects the load disturbance effects 
and provides a zero steady-state error. However, the PD plus gravity controller is not able to eliminate the load 
disturbance completely. Similar to the joint space control application results, the proposed model-free integral 


SMC controller provides robust and highly satisfactory trajectory tracking performance in task-space control 
of the PUMA 560 robotic arm. 
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Figure 11. Task-space control torque signals for each joint 
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Figure 12. Task-space tracking errors under 7; = 10 Nm load torques 


4. CONCLUSION 

A model-free continuous integral sliding mode control approach is proposed for controlling robotic 
manipulators for providing robust and high control performances. The integral term of the controller has the 
working principle of adaptive control and dynamically eliminates the effects of static terms such as gravity 
terms and load disturbances. The controller is completely model-free and utilizes the robustness feature of the 
sliding mode control for the global stabilization of robotic manipulators. The proposed controller is designed 
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for both joint-space control and task-space control implementations. The performance of the controller is 
illustrated on the PUMA 560 robotic manipulator and compared with the performance of the PD plus gravity 
control approach. The numerical results show that the model-free control method is capable of providing a 
non-overshoot response, fast settling time, zero steady-state error under load disturbance and measurement 
noise, and is robust under bounded uncertainties. 
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